Mechanism of Laser-induced Field Emission 
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We have measured electron energy distribution curves (EDCs) of the laser-induced field emission 
from a tungsten tip. Field emission from photo-excited nonequilibrium electron distributions were 
clearly observed, while no enhanced field emission due to optical electric fields appeared up to values 
of 1.3 V/nm. Thus, we experimentally confirm the emission mechanism. Simulated transient EDCs 
show that electron dynamics plays a significant role in the laser-induced field emission. The results 
should be useful to find optimal parameters for defining the temporal and spectral characteristics 
of electron pulses for many applications based on pulsed field emission. 



Applying strong electric fields (2-6 V/nm) to a metal- 
lic tip with nanometer sharpness enables field emission 
due to electron tunneling into the vacuum, producing 
continuous electron beams with high brightness and co- 
herence Illumination of such tips by femtosecond 
laser pulses has realized pulsed field emission with spatio- 
temporal control with femtosecond and nanometer reso- 
lution, making it attractive for both basic research and 
new applications like time-resolved electron microscopy, 
spectroscopy, holography, and also free electron lasers [Q- 
Despite intensive research over the last half decade 
41 6| , the emission mechanism is still controversial with 
the fundamental question being whether the optical fields 
of the laser pulse interact with the electrons in the tip as 
particles (photons) or via the electrical field proper. 

When a focused laser pulse illuminates a metallic tip, 
optical electric fields are enhanced at the tip apex due to 
plasmonic effects, and the enhanced fields induce pulsed 
field emission in combination with a moderate DC volt- 
age applied to the tip. Depending on the strength of the 
enhanced fields, different field emission mechanisms are 
considered to become dominant [1,0] • For relatively weak 
fields, single-electron excitations by single- and multi- 
photon absorption are prevalent, and photo-excited elec- 
trons are tunneling through the surface potential barrier 
{photo-field emission). On the other hand, very strong 
fields largely modify the tunneling barrier and prompt 
field emission from the Fermi level {optical field emis- 
sion). 

To reveal the emission mechanisms, measuring elec- 
tron energy distribution curves (EDCs) is the most di- 
rect method. Here we present such data, and we ex- 
perimentally confirm that the emission mechanism can 
be quantitatively described within the pure photo-field 
emission model up to optical fields of 1.3 V/nm. Simu- 
lated transient EDCs show that electron dynamics plays 
a significant role in the laser-induced field emission. 

Fig. 1(a) schematically illustrates our experimental 



setup. A tungsten tip with its axis along the [Oil] crystal 
direction is mounted inside a vacuum chamber (9 • lfV 11 
mbar). Laser pulses are generated in a Tksapphire os- 
cillator (center wave length: 800 nm; repetition rate: 76 
MHz). The temporal spread of the laser-pulse intensity 
profile is estimated to be roughly 100 fs in FWHM just 
at the tip apex. The laser light was focused to 4 /im (1/e 2 
radius) onto the tip apex. Linearly polarized laser light 
was used, with the polarization vector parallel to the tip 
axis. The tip can be heated to clean the apex and also 
negatively biased for field emission. A pinhole plate with 
phosphor coating was mounted in front of the tip to ob- 
serve emission patterns from the tip apex, and to de- 
fine a specific emission site for electon spectroscopy. A 
hemispherical analyzer (VG: CLAM2) is used to measure 
EDCs of the emitted electrons passing through the pin- 
hole. The tip can be moved along five axes as used in our 
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Fig 1: (color online) (a) a schematic diagram of the experi- 
mental setup (see text for further description), (b) a photo- 
graph of the experimental setup. The inset is the front view 
of the atomic structure of the tip apex based on a ball model. 
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Fig 2: (color online) (a) EDC of field emission from the tung- 
sten tip together with a schematic diagram of field emission 
from a Fermi-Dirac distribution. The tip voltage Vu p was 
-2300 V. (b) schematic diagrams of nonequiliburium elec- 
tron distributions just after laser absorption without (upper 
panel) and with (lower panel) electron-electron scattering, (c) 
and (d) show experimental EDCs and schematic diagrams 
of photo-field emission and photoemission, respectively. The 
magnified spectrum is also shown in (c). (c) Vu p = -2300 V 
and the laser power Pl = 50 mW. (d) Vu p = -500 V and Pl 
= 50 mW. 

previous work [j, 8]. The tip axis is set to be orthogonal 
to the pinhole plate. 

A photograph of the experimental setup is shown in 
Fig. 1(b). The field emission pattern of the clean tung- 
sten tip can be observed on the phosphor plate where 
the most intense electron emission is observed around the 
[310]-type facets. The emission sites are highlighted by 
green areas on the schematic front view of the tip apex 
in the inset of Fig. 1(b). The pinhole is positioned at 
the edge of a [310] type facet in Fig. 1(b). The position 
of the pinhole is roughly indicated by a red circle in the 
inset; the selected site is the most intense emission site 
in the laser- induced field emission 0, § . 

Fig. 2(a) shows an EDC of field emission at a tip volt- 
age Vup of -2300 V. The peak of the spectrum defines 
the Fermi energy Ep at eV. The spectrum shows a 
typical asymmetric peak, which can be understood by 
the diagram in Fig. 2(a). The field emission current 
is influenced by two factors: 1) the electron occupation 
number and 2) the transmission probability through the 
surface barrier 0, 0, H3, 13 • The occupation number is 
given by an electron distribution function f(E), which is 
the Fermi-Dirac distribution function in the case of field 



emission, and the transmission probability depends ex- 
ponentially on an area of the surface barrier indicated 
by the hatched area. Therefore, the positive energy side 
of the spectrum falls off due to a rapid decrease of the 
occupation number, while the negative energy side falls 
off because of the exponential decay of the transmission 
probability due to the increase of the surface barrier area. 
Thus a typical field emission spectrum shows such an 
asymmetric peak. An energy spread of 0.21 eV was ob- 
served, which is close to the value measured with 1 meV 
energy resolution in previous work (0.19 eV) 19], con- 
firming our reasonable energy resolution. 

In photo-field emission, the electron distribution is 
strongly modified by single-electron excitations due to 
multi-photon absorption, resulting in a nonequilibrium 
distribution characterized by a steplike profile as illus- 
trated in the upper panel of Fig. 2(b) [ll S3, E3; the 
width of each step corresponds to the photon energy hv 
(= 1.55 eV). Here we identify the step edges of one- 
photon excitation (1PE), 2PE and 3PE as shown in Fig. 
2(b). In a real situation, however, the excited electrons 
relax mainly by electron-electron (e-e) scattering on a 
time scale of a few femtoseconds, which is shorter than 
our laser pulses. As a result, the electron dynamics is re- 
flected in a smeared electron distribution, as shown in the 
lower panel of Fig. 2(b). This feature should be reflected 
in the measured EDCs. 

Fig. 2(c) shows an EDC of laser-induced field emis- 
sion at Vup = -2300 V and a laser power Pl of 50 mW. 
The spectrum shows the field emission peak undisturbed 
with identical shape and intensity as in Fig. 2(a), and 
additional peaks at the 1PE and 2PE edges are clearly 
observed. The latter show the same asymmetric shape 
as the field emission peak. Thus, photo-field emission is 
confirmed experimentally. Regarding the relative intensi- 
ties, photo-field emission from 2PE is much higher than 
that from 1PE even though the occupation number at 
1PE is higher; this is because the transmission probabil- 
ity at 2PE is quite high. Note that photo-field emission 
from 2PE has not been observed for excitation with a 
CW laser [22j]; electron distributions are supposed to be 
largely different from our case. 

An EDC at V Up = -500 V and P L = 50 mW is shown 
in Fig. 2(d). In this spectrum, the nonequilibrium 
electron distribution function becomes more noticeable. 
Field emission and photo-field emission processes are sup- 
pressed due to the low DC field, and photoemission over 
the surface barrier dominates, for which 3PE is required. 
The peak shows a spectral shape completely different 
from that of field emission. The peak maximum is lo- 
cated approximately 0.65 eV below the 3PE edge. Since 
the transmission probability is unity throughout the pho- 
toemission regime, the peak shape reflects more closely 
the electron distribution function. We observe that it is 
strongly modulated by the e-e scattering processes. 

Further investigations of the emission mechanisms 
were done by systematically measuring EDCs for vari- 
ous tip voltages and laser powers as shown in Fig. 3. All 
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Fig 3: (color online) Measured and simulated EDCs for differ- 
ent tip voltages (DC fields Fdc) and laser powers (excitation 
constant C eX c, appears in eq. (1)). All the spectra are nor- 
malized at their maximum values. The measured EDCs are 
shown by green dots with lines, and the simulated EDCs by 
red lines. Magnified spectra are also shown in the spectra at 
Vu p =-2300 V by green circles with lines for the experiments 
and red lines for the simulations. The inset (bottom left) 
shows the measured laser power dependence of intensities for 
each excitation order (Vti P =-2300 V). These intensities were 
obtained by integrating intensities over the respective energy 
regions: field emission (FE) is for -f .2 eV ~ 0.3 eV, photo- 
field emission and photoemission from fPE for 0.3 eV ~ f .85 
eV, 2PE for f .85 eV ~ 3.4 eV and 3PE for 3.4 eV ~ 4.95 eV. 
The lines are fitting curves with power functions; exponents 
are 1.3, 2.1 and 3.2 for 1PE, 2PE and 3PE, respectively. 



absorption are included in a system of Boltzmann's equa- 
tions to obtain the temporal evolution of the distribution 
function of the electron gas and the phonon gas. Fermi- 
Dirac and Bose-Einstein distributions at room temper- 
ature were assumed as initial conditions. The collision 
terms that describe e-e and e-p scattering were given in 
Refs. [2l|, |23j , and the term for the electron-gas excita- 
tion by the pump pulse was given in Ref. [24fl as 



H{E,t) = C exc F{t) B 



s {VE=fof(E-hv)[l-f(E)] 
-VE + h^f(E)[l-f(E + hu)]} (1) 



where F(t) gaU ss is a Gaussian function with a maximum 
of unity and a temporal spread of 100 fs, and C exc is 
an excitation constant which is proportional to the laser 
power and is used as a fitting parameter in our simula- 
tions. To calculate e-e and e-p scattering, we used free- 
electron and Debye models with a Fermi energy of 9.2 
eV [26[ , a Debye temperature of 400 K (2?J , and a sound 
speed for longitudinal phonons of 5220 m/s (28[. Calcu- 
lation of EDCs were done by using the Fowler-Nordheim 
theory, as in the previous works @, H, 03 ■ 

In the simulation, a laser pulse is moved in 0.2 fs steps 
from -200 fs to 1200 fs across the emission site, where 
the time zero is defined when the pulse maximum meets 
the emission site. Note that the time steps were refined 
until convergence was reached at 0.2 fs. Electron distri- 
bution functions and EDCs were calculated at each time 
step. The resulting transient EDCs were integrated over 
the entire time interval. EDCs of field emission with- 
out laser excitation were also calculated for the rest of 
one period of the laser pulse repetition cycle (approxi- 
mately 10 ns), and was added to the EDCs from the first 
1400 fs; the resulting time- integrated EDCs were normal- 
ized at the maximum intensity. Thus obtained simula- 
tions were compared with the measured EDCs normal- 
ized at their maximum intensities. Quantitative com- 



the data indicate pure photo-field emission, but domi- 
nant processes change with tip voltage and laser power. 
As a rule of thumb, with decreasing tip voltage and in- 
creasing laser power, electron emission from higher-order 
photon excitation becomes dominant. These data show a 
remarkable tunability of the emission processes via these 
two parameters. 

Throughout the whole spectra, no clear onset of optical 
field emission is observed. The inset of Fig. 3 shows a 
characteristic increase of intensities as a function of laser 
power for each excitation order at Vti P =-2300 V. The 
field emission (FE) intensity, on the other hand, remains 
constant, showing no enhanced emission due to optical 
field emission. Note that laser power values of 20, 50, 
100 mW correspond to local optical fields of 0.3, 0.9, 1.3 
V/nm on the tip apex (a field enhancement by a factor 
of 2 was included according to Refs. @, H|)- 

The influence of electron dynamics in photo-field emis- 
sion was clarified by simulating transient EDCs. The 
basic approach was the same as that used in Refs. 
[HI, Hi], |2J, [24| . e-e interaction, electron-phonon {e-p) in- 
teraction and single electron excitation by single-photon 
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Fig 4: (color online) (a) Schematic diagrams and examples of 
simulations at three time steps. Transient electron distribu- 
tion functions are shown by green lines and EDCs by red lines. 
The three EDCs are not drawn to scale, (b) time-integrated 
EDCs (upper panel) and experimental EDC (lower panel), 
(c) temporal line profile of the simulated electron pulse. 
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Table 1: Four characteristic parameters of electron pulses for 
various tip voltages and laser powers. I) Energy spread, II) 
temporal spread, III) delay of the electron pulse with respect 
to the laser pulse and IV) the number of emitted electrons per 
pulse are listed in this order at each setting. Energy spreads 
are FWHM of main photo-field emission peaks in measured 
EDCs in Fig. 3. The temporal spread and delay are taken 
from temporal profiles of electron pulses, as shown e.g. in 
Fig. 4(c). The number of electrons per pulse is calculated by 
assuming that the emission area selected by the pinhole is 10 
nm x 10 nm. 



parisons based on the absolute intensities were not done 
because the transmission function of the spectrometer is 
not well known. There are only three fitting parameters: 
1) the work function $, 2) the DC field F DC , and 3) 
the excitation constant C exc . Fitting was done first by 
adjusting these three parameters at one particular set- 
ting of the tip voltage and laser power (y tip =-500 V; 
Pl = 50 mW), and then by simulating time-integrated 
EDCs for the other settings by scaling up or down Foe 
and C exc according to the corresponding tip voltage and 
laser power. This procedure was iterated until reasonable 
fitting was obtained for all settings. 

Fig. 4(a) shows examples of the simulations at three 
time steps. The electron distribution at -100 fs shows 
a clear multiple step character, but due to the scatter- 
ing processes the distribution function becomes more and 
more smeared out as time goes on. Accordingly, the tran- 
sient EDCs change their shapes and peak positions with 
time, showing a relaxation of photo-excited electrons to 



lower energies. The resulting time-integrated EDCs are 
in good agreement with the corresponding experimental 
EDCs as shown in Fig. 4(b). From the energy-integrated 
transient EDCs one can also calculate a temporal line 
profile of the emitted electron pulse, which is shown in 
Fig. 4(c). The profile shows a delay of peak emission (31 
fs) and squeeze of temporal width (73 fs). The delay is 
caused by the electron dynamics. The squeeze is due to 
the nonlinearity of the multi-photon excitations. 

All the experimental EDCs in Fig. 3 were reproduced 
by such simulations. Fitting parameters for Foe an d 
C exc obtained from the calculations shown in Fig. 4(b) 
were scaled up or down according to the corresponding 
values of Vu p and Pl- Then the time-integrated EDCs 
were calculated and are shown in Fig. 3. Throughout 
all the various conditions, the simulations are in very 
good accordance with the experiental EDCs. Thus we 
conclude that electron dynamics play a significant role in 
photo-field emission. 

Finally, we shall comment on the relevance of these 
findings in view of designing pulsed electron sources for 
the applications mentioned earlier. From the quantita- 
tive description of experimental EDCs achieved with our 
simulations, we can be confident to predict the energy 
and temporal spreads of electron pulses, and also the de- 
lay of the emission maximum with respect to the laser 
pulses. As a rule of thumb, these values increase with 
laser power as shown in Table 1, without even consid- 
ering space charge effects in the vacuum. Such an in- 
crease is not desirable for applications in general. On the 
other hand, electron emission currents increase with laser 
power. To reach higher intensities while keeping energy 
and temporal spreads low, photo-field emission from a 
tip array [14] , |25| with an illumination by low-power laser 
pulses is considered to be the optimal solution. 
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